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Abstract. The thermal fluctuation level of the Weibel in- 
stability is recalculated. It is shown that the divergence of 
the fluctuations at long wavelengths, i.e. the Weibel infrared 
catastrophe, never occurs. At large wavelengths the thermal 
fluctuation level is terminated by the presence of even the 
smallest available stable thermal anisotropy. Weibel fields 
penetrate only one skin depth into the plasma. When excited 
inside, they cause layers of antiparallel fields of skin depth 
width and vortices which may be subject to reconnection. 
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1 Introduction 



Weibel filamentation mode by self-magnetisation of the ini- 



The Weibel fliermal ( |Weibel| \T959\ and Weibel filamenta- 
tion ( |Fried| 1 19591 instabilities, two closely related versions 
of the purely electromagnetic zero frequency instability of a 
non-magnetised plasma, have been proposed (see, e.g. 



Yoon 



and Davidson 19871 as fundamental for the generation of 



quasi-stationary magnetic fields ( Medvedev and Loebl 1 1 999| 
in high temperature or high-speed (cf., e.g., [Achterber g and| 
[Wiersma] |2007| l plasmas, respectively. The mechanism of 
this instability has been elucidated ( ,Fried 1959^ as the spon- 
taneous microscopic electron (or ion) currents flowing in a 
plasma which, under unstable conditions, become amplified 
and cause a non-vanishing magnetic field to grow in the 
plasma. As usual, the field grows up to a maximum ampli- 
tude and stabilises by some mechanism like exhaust of the 
available free energy, in the case of the thermal Weibel mode 
quasilinear depletion of the temperature anisotropy (recent 



particularly complete discussions are given in Pokhotelov et 



|aLj|2010||Pokhotelov and Amariutei||201 1^ 1, in the case of the 
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tially non-magnetised streaming particles (cf., e.g., Achter- 
|berg and Wiersmaj |2007) . The generated fields are (quasi)- 
stationary on temporal scales longer than growth and stabil- 
isation times. Nonlinear effects set on when the waves reach 
large amplitudes after amplification. At this stage, quasilin- 
ear modification of the distribution function limits further 



growth ( Pokhotelov and Amariutei 20 1 1 ) and wave-wave in 



teraction causes further nonlinear evolution (Pokhotelov etl 



al. 2010 1 which distributes particle energy into other modes. 



further limiting the growth of magnetic fields. Nevertheless, 
under non-dynamo conditions, the Weibel mechanism is thus 
an alternative to the celebrated turbulent magnetic dynamo 
mechanism, the latter being believed to be responsible for 
the generation of magnetic fields in p lanets ( ^ Christensen et. 



aL|2010 Hulotetal. 201 l| l, stars (e.g. Jones, Thompson and 
ITobias , 2010\ galaxies and on the large scale in the universe. 
The Weibel instability may grow in the universe from ther- 
mal fluctuations at early times ( [Widrow et al.j |201 l| l either 
in the chromodynamic state or later in the plasma state, and 



also after recombination during structure formation ( Ryu et 
[aT] |201 1[ ). It is believed to dominate all other mechanisms 
in magnetic field generation in relativistic shocks (cf , e.g., 
|Bykov and Treumann 201 1 ). More recently, it has been sug- 
gested to play a decisive role in the ignition of magnetic 
reconnection in collisionless current sheets (Baumjohann et 



al. 2010 Treumann et al. 2010 1. These facts, in particular 



the astrophysical application to magnetic field generation in 
shocks, let the Weibel instability attract considerable inter- 
est. Nevertheless, several open questions remain in dealing 
with it. These refer to the instability itself, independent of its 
free energy source, i.e. whether thermal or flow driven, and 
to its thermal level. In the following we very briefly clarify 
two of these: the nature of the Weibel mode, and avoiding the 
infrared catastrophe of Weibel thermal fluctuations. 
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Fig. 1. Dependence of the normalised (to maximum growth y,,,) 
Weibel growth rate 7 on wavenumber. The domain of growing wave 
numbers is limited by the skin-dpeth wave number ke = cOe/c. 



2 Weibel mode 

In which mode does the Weibel instability propagate? This 
question is non-trivial because, in unmagnetised plasma of 
density A^, electromagnetic wave propagation becomes pos- 
sible only at frequencies above the plasma frequency (Oe = 
\/ e^N I SjQnie, a condition which immediately follows from 
the unmagnetised plasma wave dispersion relation = 
]^c^ + 0)^, with k wavenumber, and c velocity of light. Zero 
frequency modes cannot propagate in this case. As a by- 
product this implies that no magnetic fields can be imposed 
into the collisionless plasma. They must have been present 
from the very beginning, i.e. they must have been born with 
the plasma. Otherwise magnetic fields can enter only a dis- 
tance of the order of the electron skin (inertial) depth, = 
c/cOg in the collisionless case. Such waves in the plasma are 
evanescent. As is well known, the situation changes, when 
the plasma is already magnetic, containing a constant mag- 
netic field B, in which case Alfven, whistler, and magne- 
tosonic waves can propagate down to extremely low fre- 
quency. Hence, what is the reason for the Weibel instabiUty 
to survive and cause real magnetic fields in plasma? 

In order to understand this point, we briefly repeat the ther- 
mal Weibel theory giving it a simple representation. Thermal 
Weibel instability requires that the plasma temperature is - 
for some unspecified reason - anisotropic, allowing to de- 
fine two temperatures 7|| , in orthogonal directions. For in- 
stance, if the plasma streams with velocity V, then the plasma 
pressure might become anisotropic. 



nkr. 



r,xi + (7;||-r,x)vv/y^ 



(1) 



(with I = 5,;, = 1,2,3 the unit tensor) thus defining 
anisotropic temperatures parallel and perpendicular to the 
flow, yielding a bi-Maxwellian equilibrium distribution. The 
free energy stored in the anisotropy is available for driving 



instability. The linear dispersion relation (see, e.g. Baumjo- 



inferred, becomes 



D(t(j,k) = (n''-ex)"e|| =0, n'' = ytV/(a^ (2) 

where ex(fi),k),e|| (a),k) are the transverse and longitudinal 
response functions. The longitudinal function is of no interest 
as it does not yield electromagnetic waves. The transverse 
response function is, using the bi-Maxwellian, given by 



ej_(co,k) = l- 



©2 1 



1- A+1 



1 



-Cz(C)]} 



0)2 



(3) 



Here Z(i^) is the plasma dispersion function of argument C, = 
(i>lk\yg^, and A — Tg^jTeL — 1 > the thermal anisotropy, 
— ITe/nie the thermal velocity, and k= (^sin0,O,A;cos0) 
is the wave vector having only two components. At very low 
frequencies o « the unstable solution obtained by putting 
a)(k) = //(k) purely imaginary, yields for the Weibel growth 
rate 



2 A:xVgx 
K koXe 



1 

yfc2 
"0 



(A + 1)(M«)^' (4) 



With A > 0, the wave will grow as long as k± < ko, where 
koXf, = VA limits the unstable range in wavenumber Maxi- 
mum growth is obtained at k,„ =ko/ y/i, close to the insta- 
bility cut-off. This is shown in Fig. [T] with maximum growth 
rate 

1(^1 



Ym (^m) " 



(5) 



thus justifying the initial assumption of small frequency, but 
the unstable waves being of wavelength somewhat larger 
than the plasma skin depth. We note in passing that the 
Weibel filamentation growth rate for two equal density 
counter streaming electron flows of speed V), becomes 



YFI=kbeVb< 1 



3e_ 
k^C^ 



2Ne 



me 
nii 



(6) 



where the index b refers to the two streaming plasmas. In the 
long-wavelength range, if the background density is large, 
Ne >Nb/2, this expression simplifies to 



Yfi 



:kXh 



(7) 



hann and Treumann 19961, from which instability can be 



where kbe — (Obe/c = A^'. This growth rate is small, unless 
Vb is large. The plot of the growth rate is given in Fig. |2] It 
shows, in contrast to the thermal Weibel mode, that the linear 
growth rate is not restricted in the first place but formally 
extends to the limit ypj /kbeVb = kXb for A: — > 00. Of course, a 
natural upper bound of k<ko is provided by the Debye wave 
number ko- 

From this one concludes that, in a strictly non-magnetised 
plasma, the Weibel mode, whether thermal or flow driven, 
is not a bulk plasma mode. Under unstable conditions when 
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Fig. 2. The growth rate of the Weibel filamentation instability as 
function of wavenumber kXj, = kX^i, ^/Nh/lNe- 



it is excited in the plasma at some location x, its magnetic 
field will be felt only in an interval [x — Af.,x + Aj,] around x, 
producing magnetic fields in a region of the size of ^2Xe- 
Therefore, in an extended plasma containing temperature 
anisotropics, the Weibel mode will necessarily cause fila- 
mentary magnetic field structures, even when excited by ther- 
mal anisotropics. These structures, because of the zero di- 
vergence requirement are located adjacently carrying Weibel 
magnetic fields of opposite direction on transverse scales of 
the electron skin depth. Such structures resemble current lay- 
ers, may be unstable against reconnection, form magnetic is- 
lands and dissipate part of the free energy by transforming 
it into plasma heat, depletion of temperature anisotropy or 
kinetic flow energy, and related effects. 

3 Thermal fluctuation effects 

The Weibel instability is not really of zero-frequency. An es- 
timate of the finite though even smaller real frequency co, 
is obtained from the real part of the dispersion relation (for- 
mally assuming cOr) under the simplifying assumption 
that CO ^ COe- This yields 



cOi ~k±c 



lA 

2~ 



-nie/mi 



A + 1 



(8) 



Since Vf,/c<C 1 this shows that the frequency is very small, 
the index of refraction is n > 1, and the wave, if propagating 
at all, must propagate on one of the low frequency magne- 
tised plasma branches, Alfven or whistler. The long wave- 
length whistler dispersion relation is 



2^2 



COCOce 



(9) 



with cOce =eB/me the electron cyclotron frequency. In order 
for the Weibel mode to propagate at such low frequency as 
a whistler (or ultimately electron Alfven) wave, this requires 
that the plasma must be weakly magnetised initially and 



(10) 



determines the wave number, with Vac the electron Alfven 
speed. This becomes possible once the Weibel instability has 
grown in the above spatial slot [jc — ,x + A,.] . 



However, initially, the plasma is not completely free of 
magnetic fields. Thermal fluctuations provide a weak mag- 
netic fluctuation background level. This has been known for 
long time ( [Landau and Lifshitzj |1959[ |1960| l for any ther- 
mal electrodynamically active medium. The explicit theory 
of thermal fluctuations in plasmas can be found in ISitenko 



( 1967 1 and Akhiezer et al. ( 1975). Application to magnetic 



fluctuations has been given by Yoon ( 2007 1 and Treumann 



|et al.| ( |2010| l. Being a thermal dielectric response property, 
thermal fluctuations can exist in regions, where real waves 
do not propagate. Otherwise, they provide a pool of initial 
rms wave amplitudes, b = y^(/?y, from which the instability 
chooses to excite a particular range of wave modes which are 
allowed to propagate. For the Weibel mode this property has 
already been exploited j Treumann et al.|[2010| l. 

On the other hand, the weak magnetic fluctuation level, b, 
also provides an average finite weak background magnetic 
field which can be chosen by the unstable wave to propagate. 
Taking this fluctuation wave field, the background Alfven ve- 
locity becomes (Va) =b/y/JlQmJVe, which enters the above 
whistler dispersion relation and fixes the wave number of the 
unstable modes which then may, indeed, propagate. 

Since the frequency becomes very low, implying that also 
k± 0, when interested in the generation of quasi-stationary 
magnetic fields, the magnetic fluctuation spectrum must be 
calculated slightly more precisely than done in the above 
references. There the magnetic fluctuation level depended 
on wave number as (fe^) °^ {kj_Xe)^^, leading to an 'infrared 
catastrophe' for vanishing k± and quasi-stationary fields. 

Such a catastrophe cannot occur in nature. A more pre- 
cise calculation should take into account two facts. First, the 
plasma also contains ions (which have been included in the 
above transverse response functions but have not yet been 
used). Second, calculation of thermal fluctuations in the ther- 
modynamic equilibrium (or quasi-equilibrium) state for the 
stable thermal Weibel mode implies 



A = 7||/rx-l <0 
The transverse response function can be written as 



e± = 1 + ■ 



ml 

«2 



1- 



(02 



(11) 



(12) 



where z — C/ V2, and <I>(z) w 2z^ for z ^ 1 . Then, in the very 
low frequency limit, O) 0, one obtains the correct ther- 
mal spectral energy density of magnetic fluctuations in the 
anisotropic-pressure plasma as 



(bHk))- 



meC^{A + \fkXe 



COe V mec^ (A+2)[fe2A2-A-m,/m,']' 



(13) 



For A > we had k^X^ < A, avoiding catastrophe and pro- 
viding a thermal background to chose for the Weibel insta- 
bility. For A < 0, the stable case, infrared catastrophe is in- 
hibited by the presence of ions until mg/m, < —A and, with 
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A: — >• 0, for —2 < A. However, when A < —2, the nature of the 
fluctuations changes. In this case magnetic fluctuations be- 
long to the firehose range which excludes the Weibel mode 
anyway. 

The above spectral energy density is a function of 
wavenumber k. It can, in principle, be used in order to deter- 
mine the wavelength range in which the Weibel mode prop- 
agates on the lowest frequency whistler (or Alfven) branch. 
This requires calculation of the rms thermal fluctuation am- 
pUtude and inserting it into the Alfven velocity. We do not 
perform this simple step here. 



4 Conclusions 

The intent of this note, which we believe having achieved, 
was a critical examination of the role of the Weibel instabil- 
ity. The Weibel instability, though known for half a century, 
has only recently been celebrated as a new and apparently 
important non-dynamo generator of magnetic fields in high 
temperature plasmas. By theory, it would provide relatively 
long-scale, 0<k±Xe < 1, quasi- stationary magnetic fields. 
The linear growth rate, however, maximises in the short-scale 
range, close to the electron skin depth (or ion skin depth in 
the ion equivalent of the Weibel mode, if excited, for instance 
in high speed ultra-relativistic flows). 

We have suggested that another restriction limits the wave- 
length, i.e. the inhibition of electromagnetic modes to prop- 
agate in dense non-magnetised plasmas. On the other hand, 
thermal fluctuations in such plasmas allow for a low level 
of rms magnetic field amplitudes which serve as background 
magnetisation. These will always be present, thereby open- 
ing a channel for the Weibel mode to propagate on the lowest 
frequency Alfven or whistler mode branches. 

At these (quasi- stationary) frequencies, thermal fluctua- 
tions were believed to undergo an infrared catastrophe. This 
is not the case. Taking into account the ion contribution 
and respecting the further condition that thermal fluctuations 
evolve in the range of stable plasma states, fluctuation the- 
ory is vahd up to infinitely long wavelengths. The magnetic 
spectral energy density behaves analytically. 

Though the weak rms magnetic fluctuation background in- 
deed solves the problem of Weibel penetration into a plasma, 
the Weibel magnetic fields are on scales not much longer than 
the skin depth. Divergence conditions then require that these 
fields form vortices of neighbouring antiparallel fields, re- 
sembling current layers. These may undergo reconnection, 
forming islands and thus limit the Weibel amplitude growth 
in various ways: by structuring the field, plasma jetting, and 
heating the plasma until depletion of the anisotropic temper- 
ature, which is the free energy source of the instabiUty. 
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